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\  ABSTRACT 

\ 

The  lower  hybrid  drift  instability  was  studied  with  a  two  dimensional  elec¬ 
trostatic  simulation  code.  Simulations  showed  good  agreement  of  the  measured 
local  growth  rates  and  frequencies  with  the  results  of  local  theory  during  the 
early  stage  of  wave  growth.  At  later  times  nonlocal  effects  become  important, 
and  a  coherent  mode  structure  develops.  This  normal  mode  was  observed  to 
propagate  up  the  density  gradient. 

At  zero  plasma  beta  and  zero  electron  temperature,  we  found  that  the 
lower  hybrid  drift  instability  is  stabilized  by  the  local  current  relaxation  due  to 
both  ion  quasilinear  diffusion  and  electron  £xf?  trapping  which  causes  electron 


heating  to  occur. 


I.  INTRODUCTION 


In  the  past  several  years  the  lower  hybrid  drift  instability1  has  attracted  considerable 
interest  within  the  plasma  physics  community,  since  it  is  likely  that  this  mode  limits  plasma 
confinement  in  theta  pinches  and  field-  reversed  configurations.  This  instability  was  simulated 
by  two  of  the  authors  previously  in  the  low  drift  velocity  regime  (  vd«  v„  )  by  using  an  one¬ 
dimensional,  electrostatic  particle  hybrid  code2 .  v  v^-  vdl  is  the  difference  between  the  elec¬ 
tron  and  ion  cross-field  velocities,  and  vrt— ( r,/m,)'A  is  the  ion  thermal  velocity.  It  was  found 
that  if  the  relative  electron-ion  drift  velocity  is  kept  constant  in  time,  which  models  a  finite  beta 
plasma3  ,  ion  trapping  causes  saturation  of  the  instability.  If  this  drift  is  allowed  to  decrease 
consistent  with  momentum  balance,  which  is  related  to  the  zero  beta  plasma  case,  then  satura¬ 
tion  is  due  to  current  relaxation4,5  .  Two  dimensional  particle  simulations  of  this  instability 
with  a  finite  plasma  beta  value  for  large  drift  velocities  (  vd»  v„  )  were  performed  by  Winske 
and  Liewer6  .  Ion  trapping  was  also  the  saturation  mechanism  in  their  simulations.  Our  object 
in  this  paper  is  to  study  the  lower  hybrid  drift  instability  at  zero  plasma  beta  in  the  low  drift 
velocity  regime  with  both  a  two  dimensional  electrostatic  simulation  code  and  a  nonlocal  theory. 

A  slab  configuration  is  used  with  a  density  gradient  in  x  (Fig.  1).  In  the  initial  Vlasov 
equilibrium,  the  ion  pressure  gradient  balances  the  zeroth  order  ambipolar  electric  force  on  the 
ions.  Since  the  characteristic  frequency  of  the  lower  hybrid  drift  instability  is  much  greater  than 
the  ion  cyclotron  frequency,  the  ions  may  be  treated  as  unmagnetized  particles. 

Simulations  show  good  agreement  of  the  measured  local  growth  rates  and  frequencies 
with  the  results  of  local  theory4  during  the  early  stage  of  wave  growth  when  the  wave  energy  is 
mostly  localized  at  the  region  xq  where  electrons  have  the  largest  £x  B  drift  velocity.  After  a 
transit  time,  i.e.,  the  time  for  the  wave  pocket  traveling  at  the  group  velocity  to  across  the 
whole  system  along  the  density  gradient.  Nonlocal  effects  become  important  and  a  coherent 
mode  structure  develops.  This  normal  mode  is  observed  to  propagate  from  xo  across  the  zeroth 
order  density  gradient  to  regions  where  the  electron  drift  velocity  vE  equals  the  wave  phase 
velocity  W Ay)  ,  and  to  be  damped  by  these  resonant  electrons.  vE  is  electron  £x  B  drift  velo- 
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city  due  to  the  ambipolar  electric  field. 

At  zero  plasma  beta  and  zero  electron  temperature,  we  found  that  the  lower  hybrid  drift 
instability  is  stabilized  by  current  relaxation  due  to  both  ion  quasilinear  diffusion4  5  and  electron 
fx5  trapping7 .  According  to  Drake  and  Huba7 ,  electrons  with  a  drift  velocity  (including  both 
the  zeroth  order  and  the  perturbed  £xfl  drift  velocity  in  the  y  direction)  which  is  greater  than 
or  equal  to  the  wave  phase  velocity,  are  in  resonance  with  the  wave,  and  can  become  trapped. 
This  electron  £x  B  trapping  was  observed  in  the  simulation.  The  trapping  caused  local  current 
relaxation  by  modificating  of  the  electron  density  profile.  The  neighborhood  of  the  point  where 
the  electrons  have  the  greatest  relative  drift  velocities  is  the  most  unstable  region  according  to 
local  theory.  In  this  region,  we  found  that  the  current  relaxes  by  modification  of  the  ion  velo- 
ci.j  distribution  function  as  well. 

In  Sec.  II,  a  description  of  the  simulation  model  and  initial  equilibrium  is  presented. 
Comparisons  of  observed  linear,  local  properties  of  the  lower  hybrid  drift  instability  with  linear 
local  theory  are  made  and  given  in  Sec.  IIIA.  Section  IIIB  is  devoted  to  nonlocal  effects  of  the 
lower  hybrid  drift  instability.  A  nonlocal  theory  is  presented.  Current  relaxation  caused  by  ion 
quasilinear  diffusion  and  electron  £x  B  trapping  is  discussed  in  Sec.  IV.  Conclusions  are  given 
in  Sec.  V. 

II.  SIMULATION  MODEL  AND  INITIAL  EQUILIBRIUM 

In  our  simulations,  a  slab  configuration  was  assumed  and  the  Vlasov  equilibria  are  func¬ 
tions  of  x  only.  All  the  zeroth  order  drift  velocities  are  in  the  y  direction  and  the  uniform  mag¬ 
netic  field  is  in  the  z  direction.  The  simulations  were  carried  out  in  the  x-y  plane  by  using  the 
two-dimensional  electrostatic,  fully  nonlinear  particle  code,  EZOHAR89  .  In  EZOHAR,  the 
boundary  conditions  along  x  are  inversion  symmetry9  at  the  high  density  side  and  a  reflecting 
boundary  at  the  low  density  side.  The  simulation  system  is  periodic  in  the  y  direction. 

Since  the  mode  frequency  of  the  lower  hybrid  drift  instability  is  much  higher  than  the  ion 
cyclotron  frequency,  ions  are  assumed  to  be  unmagnetized.  In  the  equilibrium,  it  was  assumed 
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that  the  ion  pressure  force  is  salanced  oy  the  zeroth  order  ambipolar  electric  force.  From  the 
Vlasov  equation,  the  ion  distribution  function  is  then  only  a  function  of  energy  H,  where 

H,  -  m,  v:/2  +  e<b(x)  .  (1) 

Let  the  ion  distribution  function  f,(H)  be  exponential,  as 

/,(//,)  —  C,exp[(-/n,  v:/2+e(£(jr))/  r,l  .  (2) 

Rewrite  Eq.  (2)  as  a  product, 

/,(//,)-  n,(x)g,(  v)  .  (3) 

Here 


gj(v)  -  C,exp(-m, v2/^ 7",) 

is  the  Maxwellian  distribution.  The  ion  density  profile  is 


n,(x)  -  flgexp 


—  e- 


<b(x )— <6(0) 


T, 


(4) 


(5) 


n0  is  the  ion  density  at  .v—0. 

Similarly,  the  electron  equilibrium  distribution  function  is  a  function  of  two  electron 
invariants:  energy, 

Ht  —  mt  v2/2—  e^(jr)  (6) 

and  the  guiding  center  position 


X-x-vy/o)ce  .  (7) 

Let  the  electron  distribution  function  ft{He,X)  be  given  as 

ft(He,X)  -  F(X)9xp(-m'V2/2  +  e<t>)/Tt  .  (8) 

For  a  small  electron  gyroradius  and  a  slowly  varying  ambipolar  field,  i.e.,  a, {d4>/ dx)<b~x«  1  , 

the  electric  potential  can  be  expanded  around  the  guiding  center  position  as 

<b(x)  -<*(*) +  A;Ar+ .  (9) 

ax 

where 


lx  -  x-X  -  vji o„  (10) 

is  the  electron  displacement  from  its  guiding  center.  Therefore.  Eq.  (9)  becomes 


<t>(x)  -<*(*)-  vyE(X)/<j>ct 


(11) 


Substitute  Eq.  (11)  into  Eq.  (8)  and  rewrite  Eq.  (8)  as 


/,  ( H„  X)  -  F(X)exp 


mtvl(X)/2+e<t>(X)  me 


Te 


~  Jf-  vx+  vy~  A’) J  j 


- Ne(X)gl(v ) 

where  S,(X)  is  the  electron  guiding  center  density  profile  and 


g,(  v)  -  C,exp 


is  the  electron  drifting  Maxwellian  distribution. 
From  the  Poisson  equation, 
d£(x ) 


dx 


—  tirelrtjix)  —  ne(x)]  , 


(12) 


(13) 


(14) 


and  Eq.  (5),  n,(x)  and  ne(x)  can  be  determined  by  choosing  an  appropriate  E(x)  ,  where 
n,(x)  is  the  electron  particle  density.  The  E  field  with  the  form 


E{x)  -  -2£0tanh(jr/I)sech(x/I)  ,  (15) 

which  gives  a  peak  value  as  ~E0  at  x-x0-Z.sinh~l(l)  ,  was  chosen  in  our  simulations.  There¬ 
fore,  the  ion  density  profile  is 


rtj(x)  —  noexp 


2  eE0L 


tanh  2(x/  L) 


(16) 


The  electron  density  n,  can  be  obtained  by  substituting  Eqs.  (15)  and  (16)  into  Eq.  (14). 
According  to  Eq.  (12),  electrons  have  to  be  loaded  by  following  the  electron  guiding  center 
density  profile  Nt(X)  which  can  approximately  be  expressed  in  terms  of  nt  as 


N,(X)  ®  «,(*>+ -T— T-l,  ,  <17> 

4  dx2 

for  ae~-/ n.«\  .  Hence,  the  electron  guiding  center  density  Nr(X)  is  determined  once 
ax 

E(x )  has  been  specified.  Equilibrium  profiles  of  the  normalized  ion  density  nt(x)l  no  and  the 
relative  electron-ion  drift  velocity  vE(x)  are  shown  in  Figs.  2(a)  and  (b). 

There  was  a  64x64  spatial  grid  in  the  simulated  system.  The  dimension  in  x  was  42.43xD( 


and  that  in  y  was  44.43Xo,  ,  where  Xq,  is  the  ion  Debye  length.  Time  step  0.2  was 
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chosen  and  32768  particles  were  used  for  each  species.  The  mass  ratio  mj mt  in  the  simulation 
was  100  and  1  ,  where  atp,  and  <nc,  are  electron  plasma  and  cyclotron  frequencies, 

respectively.  Ty.ptcally  vEJ  v„  was  varied  from  0.6  to  3.6,  where  vE—eEJB  is  the  maximum 
electron  £xg  drift  velocity  located  at  x-.r0 

III.  SIMULATION  RESULTS  AT  LINEAR  STAGE 

A.  Local  Effects 

Figure  3  shows  a  history  plot  of  the  simulation  electrostatic  energy  of  a  single  Fourier 
mode  at  xo  ,  where  electrons  have  the  largest  £x  B  drift  velocities.  According  to  local  theory7 , 
the  region  near  xo  is  the  most  unstable  area.  The  theoretical  local  growth  rate  is  drawn  as  a 
straight  line  .  Simulation  local  growth  rates  and  frequencies  for  different  modes  were  measured 
and  compared  with  the  theoretical  results  in  the  largest  drift  velocity  region,  as  shown  in  Figs. 
4(a)-(d)  for  vEJ  v,* 0.6  ,  0.9,  1.2  and  3.6,  respectively.  The  figure  shows  good  agreement  of 
simulations  with  local  theory  during  the  early  stage  of  wave  growth. 

B.  Nonlocal  Effects 

From  the  simulations,  it  was  found  that  nonlocal  effects  become  important  at  later  times, 
and  a  coherent  mode  structure  develops  along  x  as  shown  in  Fig.  5.  Figure  6  shows  a  snapshot 
of  a  potential  contour  plot  from  a  single  mode  simulation,  i.e.,  except  for  the  ky— 0  mode,  only 
one  Fourier  component  in  kv  of  electric  potential  was  used  to  push  the  particles.  The  ky— 0 
mode  is  necessary  in  order  to  allow  an  ambipolar  electric  field  due  to  the  charge  separation. 
Note  that  no  Fourier  transform  of  electric  potential  was  made  in  the  x  direction.  Therefore,  all 
possible  kx  modes  existed  in  our  simulations.  From  Fig.  6,  we  found  that  plasma  system  can 
be  divided  into  three  regions.  At  the  most  unstable  region  where  electrons  have  the  largest 
£x  B  drift  velocities,  the  wave  vector  k  was  observed  in  the  y  direction  as  predicted  by  local 
theory.  The  plasma  has  larger  density  in  the  region  left  of  the  maximum  £x  B  drift  velocity 
point.  It  was  observed  that  the  normal  mode  of  the  lower  hybrid  drift  instability  propagates 
toward  higher  density  with  a  mean  wave  vector  k=s  kv(—7x+7y)  that  is,  with  a  kx=—ky  .  In  the 
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second  region,  the  lower  hybrid  drift  wave  was  found  to  travel  down  the  density  gradient. 
These  phenomena  may  be  explained  by  the  following  argument. 

We  have  found  the  governing  eigenmode  equation  for  the  lower  hybrid  drift  instability  in 
a  slab  geometry  to  be 10 


*Z*k,lx)  dftt(x)/ax**k,(x)  [  2  i 


dx2  l+tote/ai"  Ne(x)  dx  ky+  ko,(x)  1  +<ojj<o}e  at-kyvE(x) 


[7 _ i_ 

'V  2  4 


-<t>k  (x)  -  0 


^  1  ~  i  ^  .  ip  I,  \A/  —  U  VIO/ 

V  2  4ir\f),(x)  l+uijJiDn  \k\v„  * 

for  vE(X)«  v„  .  <f>k  (x)  is  a  Fourier  component  in  ky  of  the  perturbed  electrostatic  potential. 


<£*^(x)  —  i/»(x)exp[— J a(x)dx]  , 
and  put  this  into  Eq.  (18)  to  obtain  the  standard  form 


d</»(x) 


—  Q(ky,at,x)tl>(x)  -0  , 


where 


Q(k.„aj,x)  -  Q,(kj,,<i>,x)  +  iQ,(ky,<o,x)  , 
Qr(ky,co,x)  -  k2+  y^p-ln(l+w £/&►£)  +  ~ 


A.p,(x)  l+utlt/uil'  w-kyVE{x ) 


as*,.*.*  >-V?rtT 


“  - -  V  2  »i,(x)  ui».  ■ 

The  perturbed  electric  potential  is  related  to  ii»  through 

*(x,y,r)  -  |l+cu£(x)/«£j  i l»(x)e'Ur>”“')  .  (24) 

The  Fourier  component  <^(x)  of  the  perturbed  potential,  was  solved  for  numerically  by  using 

Eq.  (20),  and  is  shown  in  Fig.  5(b).  Comparing  the  simulation  result  Figs.  5(a)  with  5(b),  we 
see  that  the  normal  mode  structures  from  simulations  and  theory  are  very  similar. 


In  order  to  understand  the  contribution  of  resonant  electrons  and  ions,  we  now  derive  an 


,l< 


8 


equation  of  the  wave  energy  flow.  Assume  o>q  is  the  eigenfrequency  of  the  equation 

^r- <?,(*„».*)*<*)  -  0  (25) 

S.t* 

where  <uo-w,.+  i-y  and  y— 0  .  Expanding  Q(kv,a>,x)  around  a>o  gives 

QUcv,<u,x )  —  Q(kv,u)o,x)  +  ky,cjjQ,x)(cu— a>o)  ”  Qo+  iQoT~  •  (26) 

oai  at 

where  Qo“Qiky,m o,.r)  ,  .  and  to-a>o  is  replaced  by  /d/d /  .  Substituting  Eq. 

dco 

(26)  into  Eq.  (20),  we  get 

j}T-Q^-iQo~-0  •  (27) 

9jt  or 

Multipling  Eq.  (27)  with  <b '  gives 

<l>  Qq\iIi\2-  iQo>l>  **  0  .  (28) 

or  ox 

Subtracting  Eq.  (28)  from  its  complex  conjugate,  we  obtain 

Re(Qo)^  +  2Im(?okl:  +  /#-U*^^-)  +  iIm(?o(^^-«ir^)-0.  (29) 

or  ox  Ox  Ox  Or  Or 

The  physical  meanings  of  those  four  terms  in  the  above  equation  are  the  rate  of  wave 
energy  change,  the  source  or  sink  rf  wave  energy,  the  flux  of  the  wave  energy,  and  the  fre¬ 
quency  shift  due  to  growth  or  damping  of  the  wave,  respectively  The  ratio  of  the  second  term 
to  the  first  term  gives  the  growth  or  damping  rate.  Let  us  concentrate  on  the  second  term  only 
According  to  Eqs.  (21M23),  Im@0  can  be  expressed  as 


IrnQo  - 


1  1 

lim 

ykvv  f(x)  /^T  u. 

kj)/(x)  I-i -co^Ja/c, 

y— 0 

1 

'"T" 

r? 

+ 

H 

<: 

** 

1  1 

_ —  l 

■  f  y.  )  £  I,.,  _  If  ..  /y\|  X  /  Z-  ^  ' 

k[)i(x)  [-hutjf/iult 

71  K 

y  V£\ X^OJOi  ,  Ky  V£\ X)J  +  2  l  fci  V 

(30) 


The  second  term  in  Eq.  (30)  represents  the  resonant  ions  which  drive  the  lower  hybrid  drift 
instability,  a  negative  energy  wave,  as  predicted  by  local  linear  theory.  The  delta  function  in 
the  first  term  represents  resonant  electrons  whose  £xfl  drift  velocity  equals  to  the  v  phase 
velocity  of  the  wave.  The  opposite  sign  of  the  first  term  from  that  of  the  second  term  shows 
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that  those  resonant  electrons  have  stabilizing  effects  on  the  lower  hybrid  drift  instability.  Sup¬ 
pose  that  the  plasma  has  a  density  profile  similar  to  the  density  profile  shown  in  Fig.  2<ah  and 
the  electron  Ex  B  drift  velocity  as  function  of  x  is  similar  to  the  drift  velocities  in  Fig.  2(b); 
then  there  are  two  places  in  the  system  that  electron  drift  velocities  equal  the  y  phase  velocity 
of  the  wave  excited  at  the  most  unstable  region.  One  resonant  point  is  near  the  center  of  the 
plasma  (  .v=1.2  in  Fig.  2(b)  ),  and  another  is  at  the  outer  edge  of  the  plasma  (outside  of  our 
simulation  system  in  Fig.  2(b)  ).  Therefore,  when  the  lower  hybrid  drift  wave  is  excited  at  the 
largest  drift  velocity  region,  the  wave  packet  will  travel  in  x  to  these  two  electron  resonant 
points,  where  it  dumps  wave  energy  into  the  resonant  electrons. 

IV.  SATURATION  MECHANISMS 

Saturation  mechanisms  of  the  lower  hybrid  drift  instability  in  a  uniform  magnetic  field  at 
zero  plasma  beta  and  zero  electron  temperature  were  studied.  Simulation  results  show  that  the 
lower  hybrid  drift  instability  is  stabilized  by  current  relaxation  which  is  due  to  both  ion  quasil- 
inear  diffusion  and  electron  £xfl  trapping.  According  to  local  theory,  the  neighborhood  of  the 
point  where  electrons  have  the  largest  relative  drift  velocities  is  the  most  unstable  region.  In 
this  region,  we  found  that  the  current  relaxes  by  modification  of  the  ion  velocity  distribution 
function.  Figures  7(a)  and  (b)  are  ion  velocity  distribution  functions,  averaged  in  x  over  the 
whole  system,  versus  vx  and  y,,  ,  respectively,  for  vEJ  vu~ 0.6  from  a  single  mode  simulation. 
Small  modifications  of  the  ion  distribution  function  occurred  after  saturation.  The  slight 
flattening  around  a  small  negative  vx  in  Fig.  7(a)  shows  that  the  wave  propagates  in  the  nega¬ 
tive  x  direction.  Similarly,  Fig.  7(b)  shows  that  the  wave  also  propagates  in  the  positive  v 
direction.  Hence,  the  lower  hybrid  drift  mode  propagates  up  the  density  gradient.  For  a  larger 
drift  velocity,  the  lower  hybrid  drift  wave  is  more  localized  around  the  most  unstable  region. 
Figures  8(a)  and  (b)  are  ion  distribution  functions  averaged  over  x  between  x— 21Ax  and 
x— 32Ax  ,  which  is  the  most  unstable  region  in  the  system,  for  vEJvn~3.6  from  a  single  mode 
simulation,  where  Ax  is  the  grid  cell  size  and  there  are  64  grids  across  the  system  in  x.  Larger 
modifications  of  the  ion  velocity  distribution  functions  are  shown  in  this  case. 


.  r 
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According  to  Drake  and  Huba'  .  when  the  electron  drift  velocity,  which  includes  both  the 
zeroth  order  and  the  perturbed  Ex  B  drift  velocity  in  the  y  direction,  is  equal  to  or  greater  than 
the  wave  phase  velocity,  electrons  are  in  resonance  with  the  wave.  This  electron  E x  B  trapping 
was  observed  in  our  simulations  in  the  range  that  vEtj  v,<— 0.6  to  3.6  .  Figure  9  gives  a  snapshot 
of  an  electron  density  contour  plot  for  vEJ  v.—  0.9  from  a  single  mode  simulation.  The  similar¬ 
ity  of  the  electron  density  contour  and  the  electric  potential  contour  given  in  Fig.  6  shows  that 
electrons  move  along  constant  potential  contours,  and  electrons  around  x— 1.2  are  trapped  by 
the  wave.  This  electron  ExB  trapping  causes  current  relaxation  by  local  flattening  of  the  den¬ 
sity  profile  around  the  trapping  region  as  shown  in  Fig  .  10  at  x=1.2  and  Fig.  11  at  x=0.8  for 
v£,y'v.,— 0.9  and  3.6  .  The  other  flattening  of  electron  density  profiles  in  Fig.  10  at  x=3.0  and 
Fig.  1 1  at  x=2.  are  due  to  the  ion  quasiiinear  modifications  in  the  most  unstable  regions.  Fig¬ 
ures  12(a)  and  (b)  show  electron  phase  space  (  vy  versus  x  )  for  vEJ  v.,—0.9  and  3.6  .  The  ini¬ 
tial  electron  temperature  is  zero  and  the  initial  drift  velocities  are  given  by  the  solid  curves.  We 
found  that  the  averaged  electron  current  was  reduced  after  saturation  of  the  wave.  After  sub¬ 
tracting  the  guiding  center  velocity  from  the  electron  total  velocity,  electron  heating  (shown  in 
Fig.  13)  due  to  electron  £xfl  trapping  was  observed  starting  around  the  electron  resonant 
region  (  x=1.2  as  shown  in  Figs.  2(b)  and  5  )  and  gradually  spreading  to  the  whole  plasma  sys¬ 
tem. 


Finally,  the  simulated  saturation  levels  are  compared  with  saturation  levels  predicted  by 
both  ion  quasiiinear  diffusion  theory4-5  and  electron  Ex  B  trapping  .  The  local  approximation 
was  used  in  the  ion  quasiiinear  diffusion  theory  which  gives  the  saturation  level  due  to  current 
relaxation  as 

f  l  mt  I  v£ 

nT,  8  m,  1  +w ^/<o c2,  vu  '  <31) 

where  «-<8£:>/8ir  .  Assuming  that  the  most  unstable  mode  is  dominant,  Eq.  (31)  can  be 

rewritten  as 
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where  8<t>  is  the  perturbed  electrostatic  potential.  Assuming  that  the  zeroth  order  ambipolar 
electric  field  vanishes  and  setting  the  perturbed  £x£  drift  velocity  equal  to  the  wave  phase 
velocity,  the  saturation  level  of  the  most  unstable  mode  due  to  electron  £xfi  trapping  is  given 
as 
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2  m , 


(33) 


In  Fig.  14,  Eqs.  (32)  and  (33)  are  plotted  for  100  ,  1  ,  and  kx—k,  .  Simulation 

levels  are  smaller  than  both  the  theoretically  predicted  levels  by  roughly  a  factor  of  two;  this 
could  be  due  to  the  combination  of  ion  quasilinear  diffusion  and  electron  £x  B  trapping 
occurred  in  our  simulations.  From  Fig.  14,  our  simulated  data  gives 


e&<t>  __  1 


T, 


m, 


(34) 


Furthermore,  v£—  w/ kv  varies  and  even  goes  to  zero  at  the  electron  resonant  point  in  our  simu¬ 
lations,  where  vE  is  only  the  zeroth  order  £x£  drift  velocity.  This  could  cause  electron  £x5 
trapping  occurring  at  a  lower  level. 


V.  CONCLUSIONS 

Two-dimensional  electrostatic  particle  simulations  of  the  lower  hybrid  drift  instability  in 
the  low  drift  velocity  regime  have  been  presented.  Simulations  show  good  agreement  of  the 
measured  local  growth  rates  and  frequencies  with  the  results  of  local  theory  during  the  early 
stage  of  wave  growth.  At  later  times  nonlocal  effects  become  important,  and  a  coherent  mode 
structure  develops.  We  found  that  the  lower  hybrid  drift  instability  is  stabilized  by  the  local 
current  relaxation  due  to  both  ion  quasilinear  diffusion  at  the  most  unstable  region  and  electron 
£x  B  trapping  around  the  electron  resonant  region.  Electron  heating  due  to  electron  £x£ 
trapping  was  observed  at  the  electron  resonant  region. 
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Fig.  5  (a)  Simulated  and  (b)  numerical  normal  mode  structures  of  the  most  unstable  mode 

for  100  ,  w^Jw*,**  1  ,  and  v£yv„=0.9  Both  of  them  show  that  wave  ampli¬ 


tude  vanishes  around  or— 1 . 2  which  corresponding  to  the  electron  resonant  point. 


Potential  contours  at  (D^,f“180  for  mj m,“  100  ,  Wpjn>c v«  anc* 
^yX  />,— 0.707  .  The  dashed  curve  shows  the  wave  front  of  the  wave.  Note  that  the 
vector  perpendicular  to  the  wave  front  is  roughly  at  135°  for  l<x<2  meaning 
kx~-k  .  A  small  positive  kx  component  is  also  observed  for  x>5 


region  shows  that  the  wave  propagates  up  the  density  gradient. 


Electron  density  contours  for  a  single  mode  (  k, \ £>,=0.707  )  at  u>p,r=179.8  with 
v£  '  v,=0.9  .  Note  similarity  with  Fig  6  <t>  contours,  implying  that  electrons  move 
along  the  potential  contours  and  trapping  occurs  around  jr—1.2  . 


Fig.  10  Electron  density  profile  at  (  dashed  curves  1  and  after  saturation  (  solid 

curves  )  for  m,/ mr~\0Q  .  iu^,j(ii2cr= 1  .  and  v,,=0.9  .  The  flattening  around  a  =  1.2 
is  due  to  electron  £xfi  trapping,  and  the  second  flattening  between  A  =  2.5und3.0  is 
caused  by  ion  quastlinear  diffusion  in  the  most  unstable  region 


Fig  12  Electron  phase  space  I  v,  versus  x  I  for  (3)  0  9  and  (b>  eh.  respectively 

The  initial  electron  temperature  is  zero  and  the  initial  drift  velocities  arc  given  in 
solid  curves  The  averaged  electron  current  was  reduced  alter  saturation  of  tne  wave 
around  the  largest  drift  velocity  regions 


rip  i3  Hicctron  phase  space  plots  after  suhtrncing  tne  guiding  center  velocity  from  the  total 
electron  velocity,  (a)  Kinetic  velocity  ur  versus  .v  ,  and  (b)  ur  versus  y  ,  where 
77— V —  c£x  Bl  B2  and  E  is  the  total  electric  field  in  the  system.  The  spreading  m  u, 
and  uy  shows  that  electrons  were  heated  during  the  growth  oT  the  wave  around  tr.e 
electron  resonant  region  (  ,v=1.2  ),  and  then  the  heat  is  diffused  to  the  wnoie  svs- 


Fig  14  Saturated  lower  hybrid  drift  mode  perturbed  potential  as  function  of  v„  for 

mj mr~  100  ,  top/w?,—  ]  ,  and  T,/T,— 0  .  Two  saturation  mechanisms  are  compared: 
current  relaxation  (cr)  and  electron  fxfl  trapping  (  E^B  ).  Simulation  results  are 
represented  by  dots,  roughly  half  that  of  either  theory  given  separately. 


